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ABSTRACT 

Context. Most of the hadronic jet models for quasars (QSOs) and microquasars (MQs) found in literature represent 
beams of particles (e.g. protons). These particles interact with the matter in the stellar wind of the companion star in 
the system or with crossing clouds, generating 7-rays via proton-proton processes. 

Aims. Our aim is to derive the particle distribution in the jet as seen by the observer, so that proper computation of 
the 7-ray and neutrino yields can be done. 

Methods. We use relativistic invariants to obtain the transformed expressions in the case of a power-law and power-law 
with a cutoff particle distribution in the beam. We compare with previous expressions used earlier in the literature. 
Results. We show that formerly used expressions for the particle distributions in the beam as seen by the observer are 
in error, differences being strongly dependent on the viewing angle. For example, for F — 10 (r is the Lorentz factor 
of the blob) and angles larger than ~ 20°, the earlier-used calculation entails an over-prediction (order of magnitude 
or more) of the proton spectra for E > Tmc 2 , whereas it always over-predicts (two orders of magnitude) the proton 
spectrum at lower energies, disregarding the viewing angle. 

Conclusions. All the results for photon and neutrino fluxes in hadronic models in beams that have made use of the 
earlier calculation are affected. Given that correct 7-ray fluxes will be in almost any case significantly diminished in 
comparison with published results, and that the time of observations in Cherenkov facilities grows with the square of 
the flux-reduction factor in a statistically limited result, the possibility of observing hadronic beams is undermined. 
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1. Introduction 

The study of the possible high-energy radiation from QSOs 
and MQs lacks certainty in one central aspect, the parti- 
cle composition of the jets. Is the radiation emitted dom- 
inantly produced via inverse Compton (self-synchrotron 
and/or with external fields) or via proton-proton/proton- 
photon interactions leading to subsequent meson decay? 
In the latter case, hadronic jet models have been used in 
order to assess this possibility. Most of these models rep- 
resent beams of particles (e.g. protons) linearly propagat- 
ing in a direction normal to an accretion disc, interacting 
with the matter in the stellar wind of the companion star 
in the case of MQs, or with matter in clouds, in the case 
of quasars. The protons in the beams are usually assumed 
to be distributed with a power-law which may be cut at 
high-energies, as would be the case if they are accelerated 
in Fermi processes within the jet and subject to losses. It 
is then crucial to have a correct description of the proton 
distribution as seen by the observer, since the emissivity of 
7-ray and neutrino produced by charged and neutral pions 
make use of cross sections known in that frame. In this pa- 
per we present in detail the derivation of a beam particle 
distribution as seen by the observer and compare it with 
the earlier-used expression, finding significant differences. 
Consequences are discussed. 
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2. Derivation of the particle distribution 

In what follows, primed quantities refer to quantities that 
reside in the jet frame, whereas unprimed quantities are 
used for those in the observer frame. We consider the 
Lorentz transformations for energies and momenta 

E = T(E'+j3cp\), E' = T(E- (3cp\\), 

P ,i = r(pf, + pE'/c), P \ = r(p,| - pE/c), (i) 

where p» — pcos(9) and p± = psin(#) are the parallel and 
perpendicular momentum, respectively. For the relativistic 
case, mc <C p and mc 2 <C E, the transformation of the 
energy simplifies to E = DE 1 , with D = l/[r(l-£cos(6»))]. 

The elementary invariants under Lorentz transforma- 
tions are (see, e.g., Dermer & Menon 2010) the invariant 
four-volume: dxdt = dVdt, the invariant phase-space ele- 
ment: dp/E = p 2 dpdfl/E — > ededil, with the final expres- 
sion applying to photons and extremely relativistic parti- 
cles, and the invariant phase volume: dv — dpdx, where 
bold-faces stand for three-dimensional space magnitudes. 
Because the number N of particles or photons is invariant, 
one has that 

dN _ 1 dN 

IbJ ~ pi dvdpdn ( ' 
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is also invariant. For photons and extremely relativistic par- 
ticles the last invariant becomes 



1 



dN 



e 2 dVdedn ' 



(3) 



Use of the equality E 2 = c 2 p 2 +m 2 c 4 , allows to write p 2 dp = 
(pE/c 2 )dE. Use of the latter in Eq. ^ implies that the 
expression 



1 



dN 



pE dVdEdVL 



(4) 



is invariant under Lorentz transformations. Again, for pho- 
tons and extremely relativistic particles pc = E, (= e, in the 
notation above for this case) and the last invariant becomes 
expression ([3]). We define the differential number density of 
particles as n(E,SY) = dN /[dVdEdtt], i.e., the differential 
number of photons or relativistic particles with dimension- 
less energy between E and E + dE that are directed into 
differential solid angle interval dQ in the direction of some 
physical volume dV. Because of the invariant Q and the 
definition of n(E, Q) we can write that 



n(E,n) = n'(E',n') 



P 



E 



j/ J \E' 
which is valid in general, and 
E 2 



n(E, fi) = n'(E',n') 



E' 2 



(5) 



(6) 



which is valid for photons and extremely relativistic parti- 
cles. Eq. ([5]) just says that the expression for the distribu- 
tion function /(r, p) appearing in kinetic theory of gases, 
such that the product /(r,p)dpdy is the number of par- 
ticles lying in the a given volume element dV and having 
momenta in definite intervals dp, is Lorentz-invariant (see 
e.g., Eq. 10.5 of Landau & Lifshitz 1987). Note that 



/(r,p) 



dN 



dN 



dN 



(7) 



dV dp dV p 2 dp dft {pE /c 2 )dV dE dn 
and the invariance of /(r,p) is that expressed by Eq. Q. 

2.1. The extremely relativistic case 

In this case the particle density in the jet frame, if it is a 
power-law defined as 



n'(E',n') = (A/4ir)E' 



(8) 



transforms into the observer frame, by virtue of the previ- 
ous formulae ([6j and [8]), as, 

n{E,il) = (E/E') 2 n'(E',n') = {A/Ai:)D 2 {E/D)- a (9) 
= (A/4n)E- a D 2+a . 



2.2. The general case 

Using Equations ([l} and the one for E 2 , it is possible to 
find that 



E/E' = E/[T(E-j3cpcos(6))] 

1 



r (l -/3cos(6») v /l-m 2 c 4 /£; 2 



(10) 
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Fig. 1. A power-law spectrum of protons in a beam as seen 
by the observer for different values of viewing angles, and 
an example of Lorentz factor and slope (solid lines). The 
grey-dashed lines stand for the earlier derivation results for 
the same parameters (see ^2.4). 



In addition, using the Lorentz transformation for momenta 
(Eq. [l]and the fact that p^_ = psin(Q) =p' ± ), one has 

p 12 = p 2 sin 2 (f5) + r 2 (pcos(f5) - PE/c) 2 ) 

sin 2 (6») + r 2 (cos(6<) - PE/^E 2 - m 2 c 4 ) 2 j . (11) 

Using the latter expression to get p/p' and ins erting this 
into Eq. ([5]), together with Eq. |8]) and Eq. (10) we get, 

n(E,Sl) = 



(12) 



p-a-i^-a^ _ p cos(6») v / l - m 2 c 4 /E 2 )- a - 1 



sin 2 (6>) + T 2 cos(0) 



y /l-m 2 c i /E 2 



1/2 



Note that in the limit, when mc 2 E, the latter expres- 
sion reduces to {A/Air)E~ a D 2+a . To prove the latter we 

use that (sin 2 (6») + r 2 (cos((9) - fi) 2 ) 1 ' 2 = 1/D, which can 
be established by squaring both sides, and summing and 
subtracting T 2 sin 2 (6>) + T 2 (3 2 cos 2 (6>) on the left hand side, 
with the additional use of the definition of T 2 = 1/[1 — /3 2 ]. 
Note also that we are not assuming a specific jet compo- 
sition, thus it is a general result for any kind of par ticle 
distributions. Figure [l] shows the resulting expression (12) 
for a power-law spectrum of protons in a beam as seen by 
the observer, for different values of viewing angles, and an 
example of Lorentz factor and slope (solid lines). Note that 
for high energies, when E > Trac 2 , n(E, ft) is a power-law 
distribution, as the original primed one, albeit it presents 
a different behavior for lower values of E. 



2.3. A power-law with an exponential cutoff 

We shall briefly also consider the case in which the intrin- 
sic particle distribution in the beam is a power-law with an 
exponential cutoff. This would be a natural resulting con- 
sequence of particles being subject to losses on the same 
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site where they are accelerated. As such, the cutoff appears 
in the primed referenced frame, where 



n'(E', Q') = {A/^)E'- a exp {-E'/E c 



(13) 



Note that E cut should be fixed in the beam frame, where ac- 
celeration and losses are supposed to occur, but once fixed 
it is at all effects a number, and it is not subject to Lorentz 
transformations. However, the energy variable in the expo- 
nential function is. The result is 



n(E,Q)=^-x 

■p— a — 1 Tp— a 



(14) 



2\-a-l 



E~ a (l -j3cos(6)yfi - m 2 c A /E 2 ) 



sin 2 (0) -t-T 2 cos(0) 



V_E 2 



exp ( -^-(1 -Pcos(9)y/l -m 2 c A /E 2 J . 

V ^cut / 

When mc 2 <C E, the exponential cutoff exp (— E'/E cu t) 
gets transformed into exp (— E/DE cut ), with D being the 
bulk Dopplcr factor. 



2.4. Comparison with Purmohammad & Samimi (2001) 

Purmohammad and Samimi (PS01) states that the 
isotropic particle distribution transformed from the jet into 
the observer frame results in (in this subsection and to avoid 
confusion we label their expressions with a *) : 

n4£,!!) = V +1 rx (15) 



47T 



(1 - /3cos(6») v / l - m 2 c 4 /E 2 )- a 

V2 • 



y/l-m 2 c i /E 2 



sin 2 6 + T 2 cos(0) 



Purmohammad and Samimi's (2001) paper does not pro- 
vide a derivation of their Eq. (3), quoted above as ( 15 1, but 
earlier private communication with one of us was sufficient 
to reconstruct it. The steps are the following. They start 
by assuming conservation of particles in the form 



n*(E, tt)dVdVLdE = <(£', fl')dV' dfl' dE' , 
thus one has 

dV dfl' dE' 



n*(E,n)=n' it (E'(E,6),n>(E,e)) 



dV dil dE 



(16) 



(17) 



From the Lorentz-invariance of the phase-space element 
dp/E, the invariance of pdVldE can also be proved. Then, 
one can write that 



P_ 
V' 



dE' dCl' 
iLElin' 



(18) 



Using the expressions der ived above for p/pf ', and the equal- 
ity dV = TdV, into dl7]), one gets 



T(A/4Tr)E'- a 



(19) 



\Jsm 2 {9) + r 2 (cos(6») - /3E/VE 2 - m 2 ^) 2 

From here, after replacing in the numerator E' = T(E — 
(3cpcos{9)) = T(E — (3 cos(0)V E 2 — m 2 c 4 ), the expression 
( 15 ) immediately follows. 



For mall y, one can say that there is a problem hidden 
in Eq. (16), from where the rest of the derivation follows. 



The invariant phase volume is not dVdQ.dE but rather 
p 2 dVdfldp. To avoid inconsistencies, invariants (from where 
to derive transformation laws) should be written as opera- 
tions over individually invariant factors, e.g., dN over du, 
rather than from identities (which is what one obtains re- 
placing the definition of n and n! in Eq. (16)). However, 
the latter can also yield to the correct result too if care 
is exercised. In practice, then, the problem with the prior 
derivation is in the use of the equality dV = TdV. The 
latter equation implies that dV is the proper volume dVo; 
which is clearly incorrect, given the fact that the momen- 
tum, p', of the particles in that volume is not zero. Only 
for the proper volume the referred equality is true (see, e.g., 
eq. 4.6 of Landau & Lifshitz 1987). To clarify this issue, 
let us introduce, in addition of the two reference systems 
(observer and beam), another frame Kq in which the par- 
ticles with the given momentum are at rest. The proper 
volume dVo of the element occupied by the particles is de- 
fined relative to this system. The velocities of the primed 
and unprimcd systems relative to the system Kq coincide 
by definition with the velocities v and v' which these par- 
ticles have in the systems of the observer and the beam, 
respectively. Thus, we can write 



dV = dV ^l-v 2 /c 2 , dV = dVayJl- v /2 /c 
from which one can write 



dV _ v/l-v 2 /c 2 

dV* ~ /1 ,2/2 

1 — v' /c 2 



~E' 



(20) 



(21) 



If we use Eqs. (18) and (21) in Eq. ( |16[ ), the latter trans- 
forms back into the expression of the invariant quantity 
l/(pE) x dN/{dVdEdn) and Eq. ^ follows. 

The difference between expression ( [15] ) and our Eq. ( 12 ) 
is important. Diving it one by another we obtain 



n{E,Sl) 

ME~n) 



T- 2 (l - /?cos(60-\/l - m 2 c A /E 2 )- 1 



(22) 



which is a non-trivially dependent function of T, E, and 
9. However, note that it is not a function of a. In the 
extremely relativistic case; mc 2 <C E one notices a fac- 
tor D/T of difference. Put otherwise, in the stated asymp- 
totic regime, the ratio expressed by Eq. (22) reduces to 

(n(E,n)/n it (E,Q,)) = r~ 2 (l - cos(0)\/l - 1/T 2 )" 1 , which 
is only a function of T and 9. One can see that in the case 
of 9 = 0°, n(E,n)/n+(E,n) -> 2 for growing T. For other 
values of 9, this 2D function quickly fall off. We plot the 
comparison between our derivation and that by PS01 in 
Figure [2] both as a function of viewing angle 9 (for fixed 
T = 10) and bulk Lorentz factor T (for 9 = 20°). We see 
that the differences amounts from a factor of ~ 2 in the 
case of directly-pointing beams of blazars and microblazars 
(of very low 9, in our example, 0.6°) to one order of mag- 
nitude at least for the more common cases of moderate 
inclinations (e.g., for 20°; with the difference being signifi- 
cantly larger for larger inclinations). Note that as soon as 
9 is larger than ~ 5°, the proton spectrum in the beam 
increasingly undershoots PSOl's calculation; this leads to 
an analogous correction of previously calculated 7-ray and 
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Fig. 2. Comparison of the derived expression for the spectrum of protons in beams, following Equation (22). Left: as a 
function of the viewing angle (9) for T = 10. Right: as a function of the bulk Lorentz factor (T), for 9 — 20 1 . 



neutrino emission (typically by more than on e ord er of mag- 
nitude) in all cases in which the expression ( |15| was used. 
Note also that in the low energy decline PSOFs calculation 
always overestimates significantly (by about two orders of 
magnitude) the correct yield of the proton spectrum (see 
Fig. [I] at proton energies between 1 and 10 GeV). 

3. Discussion 

The proton beam model presented by PS01 gives an expres- 
sion for a power-law proton spectrum in a beam, as seen 
by the observer. We have demonstrated that this expres- 
sion is incorrect, since it is not complying with relativis- 
tic invariants. We have derived the correct transformation 
from the jet frame particle density to the observer's, which 
strictly obeys these Lorentz invariance constraints. By com- 
paring our result with PS01, we noticed that differences are 
strongly dependent on the viewing angle. For many cases, 
PS01 calculation entails a significant over-prediction (order 
of magnitude) of the proton spectra for E > Tmc 2 , whereas 
it always significantly over-predicts (two orders of magni- 
tude) the proton spectrum at lower energies, disregarding 
the viewing angle. The 7-ray luminosity is related to the 
particle density through the corresponding emissivity of the 
processes by which the particles interact. For instance, the 
charged and neutral pion decay channels is the most effec- 



tive yield of photons and neutrinos in many situations, and 
its emissivity is directly proportional to the particle density 
n(E 7 Q). A large number of publications relying on the cal- 
culation by Purmohammad & Samimi (2001) for estimating 
7-ray and neutrino emission yields in hadronic beam sce- 
narios are thus affected in the way described above. 

One particular example can be seen in the hadronic 
MQ case (Romero et al. 2003^]and their applications (e.g., 
Romero et al. 2005, Orellana et al. 2007, Reynoso et al. 
2008). These are all affected on the light of our results. MQs 
can have r ~ 2, but 9- values above 10° are typical. MQs 
can precess, as is the case of SS433, where 9 can change even 
up to 20° around an already large mean viewing angle. 

Only a small part of the 7-ray AGNs belong to the non- 
blazar class. The majority among those associated with ra- 
dio sources possess relatively large core dominance (Abdo 
et al. 2010), indicating that their viewing angle is likely 
not large. Assuming 9 < 20°, if T < 10, PS01 deviation 
is not more than a factor of ~ 6 (at E > mc 2 ). For the 
average blazar with L = 10 and 9 < 10° the difference is 
a factor of ~ 2. However, it has been argued that during 



Note that the eq. (2) in the latter paper, has a missing 



square in the denominator factor (cos(#) — f3E/ \/~{E 2 



'))■ 



This typo has been adopted in several of the papers quoting 
Romero et al (2003). 
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extreme TeV flaring events of some blazars (e.g. Mkn 501, 
PKS 2155-304), Lorentz factors greater than r ~ 50 are 
required (e.g. Begelman et al. 2008). If the viewing angle 
remains at 9 ~ l/r qu i ot ~ 1/10 ~ 5.7° during such events, 
PS01 formulation error can grow to more than an order 
of magnitude. Furthermore, a few radio galaxies with large 
viewing angles are also detected at 7-ray energies: M87 with 
viewing angle 6 <^ 40°, and Cen A where estimates for the 
viewing angle range from 15° to 80°. These radio galaxies 
have been already modeled with the PS01 derivation (e.g., 
Reynoso et al. 2010). Here, the difference between the PS01 
formula and ours can reach up to a factor 20 for Y not larger 
than 5. 

Finally, we note that models designed to predict future 
detections of large- viewing angle sources with higher sensi- 
tivity instruments than currently available will greatly over- 
predict the number of sources if the PS01 formula is used. 
Similarly, a model for an AGN which overestimates the 7- 
ray flux could make wrong associations of unidentified 7-ray 
sources. Given that 7-ray fluxes will be in many cased di- 
minished, and that the time of observations in Cherenkov 
facilities grows with the square of the flux-reduction factor, 
the likelihood of observing hadronic beams is undermined. 
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